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Conduct  theoretical  and  experimental  investigations  to  define 
the  mechanisms  and  control  parameters  of  gas  diffusion  electrodes  in 
fuel  cells.  Work  will  include,  but  not  be  limited  to 

1.  Establish  a  mathematical  model  of  the  overall  diffusion 
process  including  terms  for  contributions  of  known  mechanisms  in 
porous  media. 

2.  Develope  expressions  for  contributing  mechanisms  which 
contain  the  physical  properties  of  the  reactants  and  products,  and 
the  physical  properties  and  geometrical  parameters  of  the  electrode. 

3.  Conduct  experimental  work  to  verify  the  mathematical  model. 
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1 . 0  SUMMARY 


Previous  work  on  mass  transport  limitations  of  porous  gas  diffu¬ 
sion  electrodes  (1)  was  extended  to  consider  the  case  where  the 
penetration  of  the  electrolyte  into  a  single  cylindrical  pore  is 
greater  than  the  radius  of  the  pore  so  that  the  liquid  in  the  pore 
can  be  considered  a  semi-infinite  cylinder. 

An  elliptical  concentration  profile  was  used  in  this  work  to 
represent  the  composition  in  the  gas-liquid  interface.  The  major 
assumptions  of  the  model  are: 

1.  Single  pore  of  uniform  circular  cross  section. 

2.  Electrolyte  penetration  much  greater  than  pore  radius. 

3.  Instantaneous,  complete,  and  irreversible  reaction  at 
surfaces  wetted  by  electrolyte. 

4.  No  reaction  at  dry  surface. 

5.  Cylindrical  symmetry. 

6.  Gas- liquid  equilibrium  at  the  interface. 

The  results  of  this  work  indicate  that  the  maximum  current 
density  supportable  by  mass  transfer  is  greater  in  the  case  of  a 
semi-infinite  cylinder  than  that  deduced  for  a  semi-infinite 

annulus.  For  an  oxygen-water  system  this  model  predicts  a  maximum 

2  -6 
current  density  of  5.1  amps/cm  at  the  optimum  radius  of  2.15  x  10 

cm;  for  the  semi- infinite  annulus  a  maximum  current  density  of  1.45 

2  *.  6 

amps/cm  occuririg  at  an  optimum  pore  radius  of  1.86  x  10  cm  was 


found . 
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It  was  found  that  for  a  pore  penetration  of  one  pore  diameter, 
dissolved  gas  reaching  the  pore  walls  wetted  by  the  electrolyte 
amounted  to  99.3%  of  the  gas  crossing  the  gas-liquid  interface. 

This  result  indicates  that  a  penetration  on  only  one  pore  diameter 
is  sufficient  to  make  the  semi-infinite  cylinder  the  effective 


model . 


2 . 0  INTRODUCTION 


The  model  of  a  porous  gas  diffusion  electrode  controlled  by- 
mass  transfer  which  was  the  subject  of  a  previous  re~p6rf  (1)  has 
been  modified  to  consider  the  case  where  the  electrolyte  penetrates 
an  appreciable  distance  into  the  pore  so  that  the  cylinder  of 
liquid  electrolyte  inside  the  pore  can  be  considered  semi-infinite 
in  length. ' 

The  approach  of  considering  the  several  steps  involved  in 
the  total  electrochemical  process  to  be  in  series  has  been  main¬ 
tained.  Furthermore,  instantaneous  and  irreversible  electron 
transfer  at  the  reactive  site,  and  no  back  diffusion  of  reaction 
products  have  been  assumed. 

The  model  under  consideration  in  this  investigation  assumes 
that  the  inner  pore  walls  become  reactive  when  they  are  wetted  by 
the  electrolyte.  Essentially  all  of  the  other  assumptions  of  the 
preceeding  report  are  considered  applicable  to  the  present  model. 

The  principal  objective  of  this  work  was  to  determine  the 
conditions  under  which  the  current  density  resulting  from  the  flow 
of  oxygen  through  a  single  cylindrical  pore  into  a  semi-infinite 
cylinder  of  non-volatile  electrolyte  (which  in  other  respects  was 
assumed  to  act  like  water)  was  a  maximum.  Another  objective  was 
to  determine  the  extent  of  penetration  which  would  validate  the 
assumption  of  the  semi-infinite  reactive  cylinder. 
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3.0  MASS  TRANSFER  IN  A  POROUS  ELECTRODE 


3.1  Description  of  the  Model 

The  model  which  was  the  subject  of  this  work  can  be  described 
as  follows: 

Pure  oxygen  flows  through  a  pore  of  uniform  circular  cross- 
section  by  a  combination  of  Poiseuille  and  Knudsen  flows  (surface 
flow  is  treated  separately),  followed  by  diffusion  through  a  plane 
gas- liquid  interface  located  in  the  pore  due  to  the  penetration  of 
the  electrolyte  into  the  pore.  The  penetration  of  electrolyte  into 
the  pore  is  sufficiently  great  so  that  the  electrolyte  in  the  pore 
may  be  considered  a  semi-infinite  cylinder;  in  principle  the  pene¬ 
tration  is  much  greater  than  the  radius  of  the  pore.  The  reactive 
area  is  assumed  to  be  that  part  of  the  pore  wall  which  is  wetted  by  the 
electrolyte.  The  oxygen  diffuses  from  the  gas- liquid  interface  to 
the  pore  walls  where  an  instantaneous  and  irreversible  electro¬ 
chemical  reaction  occurs.  It  is  further  assumed  that  the  reaction 
products  do  not  affect  the  diffusion  field  of  the  dissolved  oxygen. 

3.2  Flow  of  Oxygen  in  a  Pore 

The  flow  of  gases  through  a  pore  of  small  diameter  may  be  des¬ 
cribed  by  three  different  mechanisms  each  of  which  prevails  at 
different  gas  pressures. 

1.  Poiseuille  flow  is  bulk  flow  due  to  a  pressure  gradient  and 
occurs  when  the  diameter  of  the  pore  is  very  much  larger  than  the 
mean  free  path  of  the  molecules  of  the  gas. 
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2.  Knudsen  flow,  occurs  when  the  diameter  of  the  pore  Is  smaller 
than  the  mean  free  path  of  the  molecules  of  the  gas. 

When  the  diameter  of  the  pore  is  of  a  magnitude  comparable  to 
the  mean  free  path  of  the  gas  molecules,  both  Poiseuille  and  Knudsen 
flows  must  be  considered. 

3.  Surface  flow  is  due  to  the  migration  of  adsorbed  molecules 
owing  to  a  concentration  gradient  on  the  surface.  The  contribution 

to  the  total  flow  of  this  phenomenon  is  negligible  for  a  pore  diameter 
of  the  order  of  the  mean  free  path  of  the  molecules  of  the  fluid  (9) 
and  will  be  neglected  in  the  mathematical  analysis  of  this  model. 
However,  the  same  type  of  correction  as  the  one  discussed  in  a  pre¬ 
vious  report  (d.)  would  be  applicable  for  the  case  where  the  surface 
flow  could  not  be  neglected. 

3.3  Flow  of  Oxygen  Through  the  Interface 

Ariet  and  Distefano  (4)  cited  experimental  evidence  in  support 
of  the  thesis  that  the  resistance  of  the  gas-liquid  interface  to  the 
diffusion  of  oxygen  could  be  neglected  under  the  flow  conditions  of 
these  models.  Therefore,  the  oxygen  on  the  liquid  side  of  the 
interface  was  considered  to  be  in  equilibrium  with  the  oxygen  on 
the  gas  side. 

3.4  Diffusional  Transport  in  the  Liquid 

The  transport  of  the  oxygen  molecules  from  the  gas- liquid 
interface  to  thereactive  surface  is  caused  by  the  concentration 
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gradient  of  oxygen  in  the  liquid  inside  the  pore.  Whenever  a  difference 
in  concentration  of  a  solute  exists  between  two  points  in  a  dilute 
solution  there  tends  to  be  a  motion  of  molecules  of  the  solute  from 
the  point  of  high  concentration  to  the  point  of  low  concentration. 

This  motion  can  be  expressed  mathematically  by  Fick's  first  law 

SCA 

nAz  =  •  5ab~5I  3-‘ 

where  N,  is  the  flux  of  solute  A  in  the  z  direction,  D  is  the 
Az  AB 

diffusivity  (or  the  diffusion  coefficient)  of  A  in  B,  and  is  a 
measure  of  ability  of  A  to  move  through  B.  It  is  a  function  of  the 
size  and  shape  of  the  molecules  of  solute  and  solvent,  the  tempera¬ 
ture,  and  the  physical  state  of  the  system.  CA  is  the  concentration 
of  A  in  the  solution. 

3.5  Assumptions  of  the  Mathematical  Model 

The  steady-state  equations  of  motion  describing  the  flow  of 
oxygen  from  the  gas  phase  to  the  reacting  surface  were  formulated  for 
a  single  cylindrical  pore  based  on  the  following  assumptions: 

1.  The  distance  which  the  electrolyte  penetrates  the  pore  is  so 
much  greater  than  the  radius  of  the  pore  that  the  liquid  in  the  pore 
may  be  considered  a  semi-infinite  cylinder. 

2.  The  concentration  of  dissolved  oxygen  at  the  interface  is 

a  function  of  the  radius  defined  and  discussed  in  a  previous  report  (1). 
See  page  9,  B.C.  III. 
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3.  There  is  no  reaction  on  the  pore  surfsltface  no  t  wetted  with 
electrolyte. 

4.  There  is  instantaneous,  complete,  and  bad  irreversible  reaction 
on  the  surface  wetted  with  electrolyte. 

5.  The  density  of  the  gas  and  the  diffusiis-isivity  of  the  gas  in 
the  liquid  phase  as  well  as  the  viscosity  of  theihhe.gas  are  considered 
to  be  functions  of  temperature  only. 

6.  In  the  liquid  phase,  the  oxygen  is  trstranspoxrtted  by  diffusion 

only. 

7.  The  model  posseses  cylindrical  syiaetrJhtry,  i.  e.,  no  gradients 
in  the  angular  direction. 

8.  The  gas  phase  consists  of  a  pure  gas,  In  tlais  investigation, 
the  calculations  are  based  on  pure  oxygen. 

9.  The  electrolyte  is  non-volatile.  In  oo  ^  other  xrespects  the 
properties  of  the  electrolyte  are  assumed  to  be  9  eidentical  with  those 

of  water.  This  assumption  obviates  any  counterdh-rdiffus  ion  of  electrolyte 
vapor  in  the  gas  phase. 

10.  For  the  mass  transfer  rates  of  this  11  investigation,  at  the 
interface,  the  oxygen  in  the  gas  phase  is  in  equip-quilibr  ium  with,  the 
oxygen  dissolved  in  the  liquid.  The  validity  ofli> of  this  assumption  has 
been  demonstrated  in  an  earlier  report  (1). 

11.  The  flow  of  gas  by  surface  migration  m  is  negligible .  This 
assumption  was  shown  to  be  valid  for  the  conditiliiions  o  £  this  model  by 
the  work  of  Wengrow  (9). 
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3.6  Formulation  of  the  Mathematical  Model 

3.6.1  Flow  From  Gas  Phase  to  Gas-Liquid  Interface 


Figure  1.  Schematic  Diagram  of  Physical  Model 

/--' 

!■„  The  equation  which  describes  the  flow  of  oxygen  from  the 
bulk  gas  phase  to  the  gas-liquid  interface  is  the  same  that  was 
developed  by  Ariet  and  Distefano  (1). 
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R  =  radius  of  the  pore,  cm 
R'  =  ideal  gas  constant, 

M  =  molecular  weight  of  the  gas,  gm. 

5  b  constant  defined  by  Carman  (4) 


3.6.2  Flow  From  the  Gas-Liquid  Interface  Through  the  Liquid 

\ 

Ariet  and  Distefano  showed  that  the  equation  describing  the 
motion  of  the  gas  from  the  gas- liquid  interface,  through  the  liquid. 


to  the  reactive  site  is  Laplace's  equation  in  cylindrical  coordinates. 


d2C. 


aVi  av  . 


=  0 


3.6-2 


This  equation  is  to  be  solved  subject  to  the  following  boundary  condi-* 
tions  which  can  easily  be  deduced  from  the  assumptions  of  the  model. 

B.C.  I  At  z^O  and  r=R,  CA(R,z)  =  0 
B.C.  II  At  r^O  and  z=*»,  C.(r,°o)  =  0 

r 

B.C.  Ill  At  z=0  and  r£R,.  CA(r,0)  =  C°  1-  (|)  i 

If  we  apply  the  finite  Hankel  Transformation  ((9)  to  equation  3.6-2 
we  obtain 

(&  ’  V) a*  ■ 0 


where  C. 


'J  rC(r, 


z)  J  (a.r)dr 
o  j 


3.6-4 


and  is  defined  by  the  following  equation 


J0(a.jR>  -  0 


3.6-5 
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The  solution  of  equation  3.6-3  is 
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li 
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1  0 
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3 


C  =  Ae‘ajz 
A 


+  Be^J2 


3.6-6 


If  we  apply  B.C.  II,  then  B  =  0  and 


C.  =  Ae‘aiZ 
A 


If  we  apply  B.C.  Ill  and  the  definition  of  the  finite  Hanfcel  Trans¬ 
formation,  we  have  at  z  =  0, 


3.6-7 


rC, 


l-(f)2 


J  (a.r)dr 
o'  J 


3.6-8 


3.6-9 


and  since  this  integral  is  independent  of  z,  we  obtain,  for  any  value 
of  z, 

v  C1  r  (‘-'I*2)’  w>dr 

Now,  let  us  consider  the  definition  of  Sonine's  integral  (2),  where  n£>n>-l 

j  (xt)tn+1(l-t2)m'n‘1dt 


0  m-n 

J  (x)  =  - — 

“  2m'nr(n 


3.6-10 


dr 


and  if  we  let  t  =  —,  dt  =  — ,  a^R  =  x,  n  =  0  and  m  =  3/2,  then 

K  K.  J 

equation  3.6-10  becomes 

J3/2<V>  '  23?f>3/2  /Rj°(a2R  t>  t  (»-(f  >2)^ 

2iU?OID 


2(aiR)3y'2  pR  /  r 

J3/2(“iE>  '  23/2r(3/2)R2i  V‘(* 


)  )  JQ(a^r)dr 


3.6-11 


3.6-12 
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If  we  compare  the  integrals  in  equations  3.6-9  and  3.6-12  we 
see  that  they  are  equal,  therefore 

C°e‘alZ  23/2  T  (3/2)  R2 

\  =  3/2  ^3/2^ai^^  ’ 

A  2(ajR  y'  ‘  J 

Now  we  can  apply  the  definition  of  the  inverse  finite  Hanke.l  Trans¬ 
formation  ($0 


to  equation  3.6-13  we  obtain 


-4  7-* 

R“  L->  .T_  ( 


J  (a.r) 


^  J1(ajR) 


,2  A 


C*. 


CA 

A 


00  J0(a.r)  C°  e'V  23/2  T  (3/2)  R2 


I 


J1(ajR) 


2(ctjR) 


3/2 


J3/2(ajR) 


or 


c  =  23/2  r  (3/2)  cA  >  -  /? 

A  A/L  (ajR)J/“ 


o  V  e  ajZ  Jpfajg)  J3/2  (ajR) 


j=1 


VajR) 


Since  (2)  T  (3/2)  =  (ttM)'2 
aiiv.  since  (2) 


J3/2^ajR)  Va.Rir/  V  a.,R 


2  X1  /sin  al  R  -  cos(ajR) 


then  equation. 3. 6- 16  becomes 


k  X 


-a+z  J  (a.r)  .sin  a.R 

o  i  f _ i_  _  cos  a^R 


JqCa.R) 


J 

V  ajR 


Differentiating  partially  with  respect  to  r  we  obtain 

n  ajR 


~V  CO 

OC.  r-, 

.  A  _  ,  „o  V 

-a -5  2  - 

e  J 

Jx(ajr) 

ai 

°r  A  ^ 

j=l 

(ajR)2 

J1(ajR) 

2 

a^R 


3.6-13 


3.6-14 


3.6-15 


3.6-16 


3.6-17 


3.6-18 


cos  a^R 


3.6-19 
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and,  if  we  evaluate  this  function  at  r  =  R,  we  obtain 


2  C 


cc.z  (-cCi)  /sin  cc.R  „  _ 

-  '  3  cos  ajR 


i  r=R 


^2  VV>  V  ajR 

j=i  3 


3. 


Since  the  total  flow  in  the  r  directiorr'at  r  =  R  from  z  =  0  to  z  =  z 


xs  siven  dv 


!  z 


W 


!  z=0 


Ar 


27fR  dz 


!r=R 


and  since  K. 


=  ~D, 


3. 


3. 


r=R 


i  r=K 


then 
!  z 

W  !  =2  G°27tR( 

iz=0  " 


(-a,) 


A  (o;ill)i'J1(a..R)  \  ajR 
3  =  a.  J  x  3 


W 


-v  I  r2, 

■2  v—i  / 1  “CL a  z v  /  sin  a.-R 

|  =47TRC°  h  Y  — A-a.-.J-L-/ 

f-,=0  ^  Z_j  ('/v  .t!  3“Y,  fr-y  sA  “'j“ 


-sin  C6jR  cosa  R\rz_a,  , 
-  3  )/  e  3  dz  . 


(z=0 


si  (otjR)~J1(<XjR) ' 


cos  a.R 
J 


-•  -.vu  C  D  ,,  V 


sin  a.R 
_ 3 _ 


cos  0..-R 


v; 

oAj-  (o^R) JJ1(ajR)  (a^R)  zJ1(ajR). 


"■CC  2  Z  « 

i  L  c:_n  a.K  e  J  cos  cCjk 


(otjR)  (ajR)A1(ajR)' 


) 


”e  can  see  that 
a-,  ea  which  is  t: 
infinite  dist: 


.as t  tr  s  terms  of  this  series  become  zero  when 
.  ..a. an  the  electrolyte  penetrates  the  pore  an 


CO 


W 


\ — i  “  sin  Co . k 

,  ...  \  i  . __J._ 

~V .  V.  VJ 


;  2— u 


,At(o:.R)  Jx(a,R)  (a,R)  j^a-jR) 


cos  cc.R  -i 

_ ] _ I 

2-  '-,«J 


6-20 

.6-21 

.6-22 

5.6-23 

.6-24 

.6-25 


3.6-26 
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and  since  by  equation  3.6-5,  (a^R)  are  defined  to  be  the  zeros  of 
the  Eessel  function  of  order  zero,  let  OjR  =  xn,  then 


W 


CO 


=  471RC 


lz=0 


:t  dab  X 


sin  x 


cos  xn  -| 


“Xn'Jjl(xn)  xn2jl(xn>^ 


3.6-27 


This  series  was  evaluated  ir>  the  following  manner;  the  first 
ten  terms  ware. summed  by  an  exact  computation  from  a  table  of  the 
zeros  of  JQ .  The  terms  from  n  =  11  to  n  =  20  were  summed  by  applying 
the  following  approximations  which  were  found  to  be  extremely  accurate 
for  the  purpose  of  this  investigation  (less  than  1%  error). 

Eor  n>l0(6) 

W  =  @  cos  (:<n  •  3.6-28 


and  for  n>10  ("5  ) 


“n  “  "  P 


hence  ^(s*)  =  ^  ^-1/4)  V'3  |  "  4>J 
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cos 


J-  (x  )  =  — 
1  ir  7 r 


Tr(n-l) 


(n-I/4) 


1/2 


2-  (-1)  n~1 
^  (n-l/4)1/2 


3.6- 29 

3.6- 30 

3.6- 31 


Now  cos  x  =  cos  (7T(n-T>)  =  (-l)n  cos  j  =  (-l)r 


sin  x  =  sin  7r(  n-1/4) 
n  1 


=  (-D11'1  sin  l  =  (-l)0'1  f 


then  equation  3.6-27  becomes 
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3.6-34 


14 


!  °o 


z=0 


■  4*ca  da» 


w 


z=0 


fo. 54261  +  V  r - - - 77T  +  - i - r 

C  2TT2nixi  L(n-1M)5/“  (n-l/4)3/2. 

:A  °AB  f0-54261  +  -Tj  (0-01298  +  0.56571) 
t  27T 


3.6-35 


^  Ff— 

271  4ri  Mn-l/4 


_ZL. 


2^^  L(n-l/4)5/2  (n-l/4)3/2JJ 


tU 
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'3.6-37 


For  the  determination  of  the  infinite  summation  in  equation  3.6-37, 

the  Euler-Maclaurin  expansion  is  invaluable.  This  expansion  is  given  by 

Miller  (^)  as 
z-l 


2-1  m-1 

Z «■> -i  ‘«*  *Z  ft  [ 

x=a  a  k=a  L 


Dk_1  f(z)  -  Dk-1  f (a) 


+  R(m)  3.6-38 


where  Bk  =  Bernoulli  Numbers  of  the  first  kind 

k-1  rik—  1 

D1  r  (a)  =  ^  f(x) 
dx 


x=a 


R(m)— 5-0  as  m— 


In  the  case  of  z  approaching  infinity,  with  f(x)  =  — - - ,  where 


n  is  positive. 


(x-k) 


2m- 1 


Z  ft  °k'1  £(z> 


k=a 

is  zero,  and  the  above  expansion  can  be  written 
00  00 


i  ■  r *«*  - 1  £  d*-£  ««.) 

x=a  ^  k=l 
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3.6-39 


15 


Therefore  the  first  term  of  the  infinite  summation  in  equation  3.6-37 


can  be  written  as 
oo 


V  1  _  f  1  ,  V  -k  nk'1  1 

^l(n-!/4)5/2  'il  (n-l/4)5/2  k'*  (n-l/4)5/2 


3.6-40 


n=21 


The  Bernoulli  numbers  of  the  first  kind  are  given  in  Miller  (7.)  as 
1 


B1  -  -  2 


B2  6 


B„  =  B„  =  B„  =  —  =  B„,  , ,  =  0 
B4  =  '  30 


B6  42 


etc . 


Equation  3.6-40  becomes 
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3.6-41 


The  expansion  on  the  left  hand  side  of  the  equation  results  in 
an  alternating  series  which  converges  very  rapidly.  Evaluating  equation 
3.6-41  at  n  =  21  results  in. 


>21 
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n=21  (n_1/4) 


jj2  =  0.00708  +  [± - - 
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=  0.00708  +  0.000256  =  0.00734 


Similarly, 
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3/2 


=  0.439  +  0.005  =  0.444 


3.6-45 


7T 


n4i 

hence  equation  3.6-36  becomes 


3/2 


=  1.39487 


z=0 


=  47TRC.  Djt, 
A  AB 


0.57193  + — r  (0.00734  +  1.39487) 
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If  we  rewrite  equation  3.6-25  as 

sin  (cijR) 
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“  “CA  V/ 
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it  can  be  seen  that  the  total  flow  to  the  walls  of  the  pore  is  equal 

to  the  flow  through  the  gas- liquid  interface  only  as  z  approaches 

~(X«  z 

infinity,  i.e„,  when  e  J  approaches  zero,  the  flow  through  the 

00 

gas- liquid  interface  equals  W^_q. 

Although  the  first  summation  in  equation  3.6-27  converged  so 
slowly  that  the  application  Euler-Maclaurin  expansion  was  necessary 
to  perform  the  summation,  the  second  summation  converged  very  rapidly 
because  of  the  term  e  ^ajz\  This  second  summation  was  carried  out 
term  by  term,  and  for  values  of  z  greater  than  ohh  pore  radius,  four 
terms  were  sufficient  to  represent  the  infinite  summation. 

Values  of  the  flow  that  reach  the  pore  walls  were  calculated 
for  different  values  of  z,  and  the  results  are  shown  in  Figure  2. 

It  can  be  seen  that  when  z  is  equal  to  2R  (electrolyte  penetration 
of  one  pore  diameter)  99.3%  of  the  flow  through  the  interface  reaches 
the  inside  of  the  pore  wall  watted  with  electrolyte.  This  is  very 
significant  for  porous  gas  diffusion  electrodes  for  it  means  that 
a  pore  penetration  by  the  electrolyte  of  one  pore  diameter  is 
sufficient  to  make  the  semi-infinite  cylinder  the  effective  model. 

3.7  Maximum  Current  Density 

The  flow  of  gas  through  the  gas  phase  and  through  the  gas- liquid 
interface  are  equal  because  they  are  in  series.  Since  the  area  that 
flux  is  based  on  is  taken  to  be  the  same,  viz.,  the  cross-sectioned 
area  of  the  pore,  the  flux  through  the  gas  phase  can  be  equated  to 
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J 

1  System:  Oxygen  Diffusing  Through  Water 


Figure  2.  Fraction  of  Gas  Flow  Through  Gas-Liquid  Interface 
Reaching  Pore  Walls 
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The  flow  of  ciygan  through  the  interlace  was  shown  to  be  equal 
to  the  total  flow  to  the  core  wall  in  an  infinite  distance  from  the 
interface.  The  flux  of  oxygen  through  the  interface  is  therefore 
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Algebraic  manipulation  of  equation  3.7-6  results  in 


5.58  x  10 


•n  Ji. 

i-p. 

t 


1.29  x  10J  R  (1+P  )  +  12.1 


3.7-7 


Under  the  restriction  that  P.  and  R  be  related  as  in  equation 

at 


3,7-7.  we  obtain  a  maximum  flux  to  the  pore  walls  (N.„) 

’  r  AR  max 


z=0 


an  optimum  radius  of  2.15  x  10  cm.  The  interfacial  pressure 
corresponding  to  this  radius  was  computed  to  be  0.51  atm.  The 
maximum  flux  to  the  pore  walls  under  these  conditions  is,  from 
equation  3.7-4, 


^AR^max 


=  1.32  x  10  5  moles/cm^.sec. 


z=0 


Therefore,  the  maximum  current  density  which  can  be  supported  by 
mass  transfer  to  a  reactive  site  is  (assuming  4  electrons  per  mole¬ 
cule  of  oxygen) 


^m-max  ^AR'max 


mol 


3.86  x  10 


z=0  cm  .sec. 


5  amps  sec 
mol 


=  5.1 


amps 

2 


cm. 


4.0  DISCUSSION  OF  RESULTS 


The  value  for  the  maximum  current  density  predicted  by  this 
model  is  approximately  three  times  greater  than  the  value  obtained 
by  Ariet  and  Distefano  under  the  assumption  of  no  electrolyte  pene¬ 
tration  in  the  pore  (semi-infinite  annulus)  (1).  This  higher  value 
may  be  explained  by  the  fact  that  the  liquid  phase  diffusion  path 
is  shorter  when  the  electrolyte  penetrates  into  the  pore  than  when 
the  gas  has  to  diffuse  into  the  liquid  outside  the  electrode  to  get 
to  the  annular  reactive  area.  Furthermore,  it  is  conventional  to 
calculate  the  current  density  based  on  the  geometric  area  of  the 
electrode.  The  area  considered  in  this  case  is  only  the  cross- 
sectional  area  of  the  pore  whereas  in  the  case  of  no  penetration 
the  area  is  the  cross-sectional  area  of  the  pore  plus  the  annular 
reactive  area. 

One  of  the  major  assumptions  of  this  model  is  that  the  liquid 
penetrates  the  pore  to  such  an  extent  that  the  pore  can  be  considered 
a  semi-infinite  cylinder.  An  important  question  immediately  arises: 
What  depth  of  penetration  relative  to  the  pore  radius  is  necessary 
to  constitute  a  semi- infinite  cylinder?  From  Figure  2  it  can  be 
-seen  that  at  an  electrolyte  penetration  of  one  pore  radius  92.4  per 
cent  of  the  flow  through  the  gas-liquid  interface  reacts  on  the 
pore  wall.  At  a  penetration  of  only  one  pore  diameter,  99.3  per  cent 
of  this  flow  reacts  on  the  pore  wall.  It  was  concluded  from  these 
results  that  any  penetration  depth  greater  than  one  pore  diameter 
is  sufficient  to  constitute  a  semi-infinite  cylinder. 
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The  assumption  that  the  presence  of  the  reaction  products  on 
the  reactive  sites  does  not  alter  the  diffusion  rates  of  the  oxygen 
in  the  liquid  will  produce  a  higher  value  for  the  current  density 
than  the  true  physical  model  in  which  these  reaction  products  will 
probably  be  a  serious  obstacle  to  the  diffusion  of  oxygen  in  the 
liquid  inside  the  pore.  It  should  be  obvious  that  as  the  wet-'.:  •  •• 
proofing  agent  is  oxidized  and  the  liquid  penetrates  in  a  pore 
whose  inner  walls  are  not  activated,  the  diffusion  rate  of  the  oxygen 
to  get  to  a  reactive  site  will  be  much  slower  because  the  gas  will 
have  to  diffuse  a  longer  distance  in  the  liquid  where  the  diffusion 
coefficient  is  much  lower  than  in  the  gas  phase.  This  will  result 
in  a  decrease  in  the  current  density.  This  effect  is  commonly 
referred  to  as  the  "drowning”  of  the  electrode.  When  the  inner 
walls  of  the  pore  are  activated,  the  penetration  of  the  liquid  in 
the  pore  will  produce  an  increase  in  the  current  density  according 
to  the  results  of  this  study.  This  seems  to  indicate  that  it  is 
desirable  to  activate  the  inner  pore  walls  of  porous  electrodes 
rather  than  to  activate  only  the  outer  surface  of  the  electrodes. 
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NOMENCLATURE 


A  Represents  component  A  in  the  binary  system  of  A  and  B 

B  Represents  component  B  in  the  binary  system  of  A  and  B 

B,  Bernoulli  numbers  of  the  first  kind  of  order  k 

k 

3 

Concentration  of  component  A  at  any  point  (r,z),  mol/cm. 

_  ^ 

CA  Concentration  of  component  A  at  any  point  (r,£),  mol/cm. 

o  3 

C^  Concentration  of  component  A  at  the  point  (0,0),  mol/cm. 

o 

Diffusivity  of  A  through  B,  cmT/sec. 

H  Henry's  Law  constant,  atm. 

J.  Bessel  function  of  the  first  kind  of  order  1 
1  J 

L  Length  of  pore  not  penetrated  by  the  liquid 

M  Molecular  weight  of  the  gas 

n  Defined  by  equation  3.6-29 

^Ar  Molar  flux  of  component  A  in  the  r-direction  at  any  point  (r,z), 
moles/cm. ^sec . 

N.J  Z  Total  molar  flux  of  component  A  to  the  reactive  pore  site  in 
e=0  the  region  z=0  to  z=z 

^Az  Molar  flux  of  component  A  in  the  z-direction  at  any  point  (r,z), 
moles/cm. ^sec. 

N.  .  Total  molar  flux  of  component  A  through  the  gas-liquid  interface, 
21  mol/cm. 2sec. 

^AzPG  H°^ar  flux  component  A  in  the  gas-phase  area  of  the  inside 
of  the  pore,  mol/cm. ^sec. 

ph  + 

P  Average  pressure  in  the  pore,  ^ k  ,  atm. 

Pg  Bulk  gas  pressure,  atm 

P^  Gas- liquid  inter facial  pressure,  atm 


i 
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X  Denotes  radial  direction,  cm 

R  Radius  of  the  pore,  cm 

R'  Gas  constant  in  appropriate  units 

t  Denotes  time,  sec 

T  Denotes  temperature,  °K 

W  |z^o  Total  moles  of  flow  to  the  pore  walls  in  the  region  z=0  to 
z=z,  moles/ sec. 

x^  Defined  by  equation  3.6-29 

z  Denotes  axial  direction,  cm 

SUBSCRIPTS 

A  Denotes  the  component  A  in  the  binary  mixture  of  A  and  B 

B  Denotes  the  component  B  in  the  binary  mixture  of  A  and  B 

b  Denotes,  the  bulk  gas 

G  Denotes  the  gas 

i  Denotes  the  gas- liquid  interface 

j  Denotes  the  root  of  the  first  order  Bessel  function  of 

the  first  kind,  J^(ajR)  =  0 

k  Denotes  the  k’"'1  order  Bernoulli  number  of  the  first  kind,  B, 

k 

n  Denotes  the  n1-11  root  of  the  first  order  Bessel  function  of 

the  first  kind,  J^(xn)  =  0 

P  Denotes  the  pore 

R  Denotes  the  reactive  pore  area 

r  Denotes  the  radial  direction 

z  Denotes  the  axial  direction 
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GREEK  SYMBOLS 

5  Constant  defined  by  Carman  (3) 

7)  Viscosity  of  the  gas,  poise 

Ctj  Defined  by  equation  3.6-5 

p  partial  pressure  of  oxygen,  atm 


27 


REFERENCES 


1.  Ariet,  M.  and  Distefano,  G.P.,  Summary  Report  No.  8,  University 
of  Florida,  March  15,  1963. 

2.  Bowman,  Frank,  "Introduction  to  Bessel  Functions",  p. 

Dover  Publications,  Inc.,  New  York  (1958). 

3.  Carslaw,  H.S.  and  Jaeger,  J.C.,  "Conduction  of  Heat  in 
Solids",  Oxford  Press,  London  (1959). 

4.  Carman,  P.C.,  Proc.  Roy.  Soc.,  203  A,  55,  London  (1950). 

5.  Janke,  E.  and  Emde,  F.,  "Table  of  Functions",  4th  ed.,  p. 

Dover  Publications,  Inc.,  New  York  (1945). 

6.  Mickley,  H.S.,  Sherwood,  T.K.  and  Reed,  C.E.,  "Applied 
Mathematics  in  Chemical  Engineering",  p.  176,  McGraw-Hill 
Book  Company,  New  York  (1957). 

7.  Miller,  K.S.,  "An  Introduction  to  the  Calculus  of  Finite 
Differences  and  Difference  Equations",  p.  113,  Henry  Holt 
and  Co.,  New  York  (1960). 

8.  Perry,  J.H.,  "Chemical  Engineer's  Handbook",  3rd  ed.,  McGraw- 
Hill,  New  York  (1950). 

9.  Sneddon,  I.N.,  "Fourier  Transforms",  p.  104,  McGraw-Hill  Book 
Go.,  New  York  (1951). 

10.  Wengrow,  Henry,  Summary  Report  No.  6,  University  of  Florida, 
December  15,  1962. 


ADDENDUM  I 


Computation  of  Concentration  Profiles 

Concentration  profiles  as  a  function  of  r/R  with  z/R  as  a 
parameter  are  shown  in  Figure  3.  The  values  of  C^/C°  were  cal¬ 
culated  from  a  truncated  form  of  the  infinite  series  in  equation, 

(3).  This  equation  resulted  in  a  convergent  series  containing 
alternating  terms.  In  all  cases,  six  or  seven  terms  were  sufficient 
to  insure  that  the  truncated  series  was  an  adequate  representation 
of  the  infinite  series. 
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ADDENDUM  II 


The  System  Air-Water 

It  was  shown  in  a  previous  report*  that  for  an  air-water  system, 

the  oxygen  flow  through  a  stagnant  layer  of  nitrogen  in  the  pore  can 

be  calculated  using  the  partial  pressure  differential  of  oxygen  as 

the  driving  force.  When  the  bulk  gas  pressure  of  air  is  P,  =1.0  atm 
.  b 

the  bulk  partial  pressure  of  oxygen  is  approximately  p^  =  0.2  atm. 

If  all  other  conditions  are  identical,  and  the  contribution  of  Poiseuille 
flow  is  neglected,  equation  3.7-5  becomes 

NAzPG  =  12-°5  (°-2  -  P±>  <la> 

Equating  equations  3.7-4  and  (la)  results  in 

-11  P± 

12.05  (0.2  -  p.)  =  5.58  x  10  (2a) 

Under  the  constraint  of  equation  (2a),  the  optimum  pore  radius  is 

2.15  x  10  6  cm,  with  a  corresponding  interfacial  partial  pressure 

p.  =0.1  atm.  Under  these  conditions  the  maximum  current  density 

2 

obtained  is  1.0  amp/cm  . 


*  Ariet,  M.  and  Distefano,  G.P.,  Summary  Report  No.  8,  University 
of  Florida,  March  15,  1963. 
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